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FREE- SURFACE  BOUNDARY  LAYER  AND 
THE  ORIGIN  OF  BOW  VORTICES 

I.  INTRODUCTION 

The  requirements  of  vanishing  tangential  stresses  at  a  curved  free 
surface  imply  a  nonzero  vorticity  which  is  convected  and  diffused  into  the 
fluid.  At  sufficiently  high  Reynolds  numbers,  the  vorticity  is  confined  to  a 
thin  boundary  layer  along  the  free  surface.  The  oscillatory  boundary  layer  at 
the  free  surface  of  travelling  and  stationary  water  waves  has  been  considered 
by  Longuet-Higgi ns  (1953,  1960),  particularly  to  determine  its  effect  on  the 
mass-transport  velocity  in  the  fluid  outside  the  boundary  layer.  Batchelor 
(1967)  has  also  described  the  origin  of  the  free-surface  boundary  layer  and 
presented  two  applications,  namely  the  drag  of  a  spherical  gas  bubble  rising 
through  a  liquid  and  the  attenuation  of  gravity  waves  due  to  viscous  effects 
at  the  free  surface.  A  free-surface  boundary  layer  is  also  present  ahead  of 
an  object  in  motion  due  to  the  surface  elevation  above  the  undisturbed  level 
and  the  associated  curvature. 

The  possible  importance  of  the  free-surface  boundary  layer  ahead  of  a 
body  moving  at  the  surface  of  a  liquid,  e.g.,  at  the  bow  of  a  ship,  was 
pointed  out  by  Landweber  and  Patel  (1979).  A  number  of  flow-visualization 
experiments  have  shown  the  existence  of  vortices  under  the  free  surface  at  the 
bow.  Such  vortices  were  first  reported  by  Suzuki  (1975)  and  Honji  (1976)  in 
two-dimensional  flow  ahead  of  a  semi -submerged  circular  cylinder.  Kayo  et  al . 
(1982)  have  repeated  these  experiments  and  confirmed  the  occurrence  of  these 
vortices.  In  three-dimensional  flow,  the  experiments  of  Kayo  and  Takekuma 
(1981)  and  Shahshahan  (1982)  with  towed  ship  models,  those  of  Kayo  et  al. 
(1982)  with  towed  vertical  cylinders,  and  some  ongoing  observations  made  by 
the  present  authors  with  fixed  cylinders  in  a  hydraulic  flume  show  horseshoe 
vortices  forming  ahead  of  the  bow  just  below  the  free  surface.  Although  bow 
vortices  are  also  discussed  in  several  recent  papers  (see  Maruo,  1983;  Mori, 
1984;  Takekuma  and  Eggers,  1984),  the  precise  mechanism  responsible  for  their 
formation,  and  the  role  they  play  in  the  breaking  of  bow  waves,  is  not  yet 
understood. 
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The  two-dimensional  free-surface  boundary  layer  is  the  subject  of  the 
present  paper.  The  necessary  conditions  at  a  curved  free  surface  of  a  viscous 
fluid  in  motion  are  examined  and  it  is  shown  that  surface  tension  plays  a 
critical  role  in  determining  the  real  flow  ahead  of  an  obstacle.  In 
particular,  the  boundary  conditions  and  the  equation  of  continuity  lead  to  a 
criterion  for  the  occurrence  of  a  stagnation  point  at  the  free  surface,  which 
may  be  identified  with  the  existence  of  a  vortex  further  downstream.  The 
theory  is  in  approximate  agreement  with  the  experimental  observations  in  two- 
dimensional  flow  noted  above,  and  may  explain  the  origin  of  the  bow  vortices. 

The  equations  of  the  free-surface  boundary  layer  are  then  derived  and  an 
approximate  integral  method  of  solution  is  presented.  This  leads  to  an 
estimate  of  the  momentum  thickness  of  the  boundary  layer. 


II.  NAVIER- STOKES  EQUATIONS  AND  BOUNDARY  CONDITIONS 


Consider  a  two-dimensional  obstacle  in  a  uniform  stream  of 
velocity  U^,  as  shown  in  figure  1.  The  vertical  distance  above  the 
undisturbed  level  far  upstream  is  denoted  by  y  and  the  elevation  of  the  free 
surface  above  this  level  is  c.  It  is  convenient  to  choose  a  curvilinear 
orthogonal  coordinate  system  (s,n)  in  which  s  is  along  the  free  surface  and  n 
is  normal  to  it.  If  the  corresponding  velocity  components  are  denoted  by  u 
and  v,  then  for  laminar  flow,  the  equation  of  continuity  and  the  Navi er-Stokes 
equations  in  the  (s,n)  directions  may  be  written,  respectively. 


(1) 


iy. 

+  K 12  an 


(2) 
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where  p  is  density,  v  is  kinematic  viscosity,  g  is  acceleration  due  to 
gravity,  p  is  pressure. 


h  i  -  1  +  <  n 

i  3hi 

_  1 _ 1  _  tc 

*  12  ~  h  '  9n  1  +  ten 


(4) 

(5) 


and  ic  is  the  curvature  of  the  free  surface. 

The  normal  and  tangential  stresses  in  the  fluid,  t  and 

nn 

xsn,  respectively,  are 


_  ,  o  9  V 

T  =  -  p  +  2  p  — — 

nn  K  9n 


JL  9  ,U  N  ,  1  9  V  i 

Tsn  p  an  ^h  )  +  h^  9s 


(6) 

(7) 


where  u  (=  pv)  is  the  coefficient  of  viscosity.  The  only  component  of 
vorticity  is 


=  tf?  “  In  (hlu)J 


(8) 


At  the  free  surface,  we  have 


v  =  0 
o 


(9) 


and  since  the  tangential  stress  vanishes,  equation  (7)  gives 


/9u^ 

'9n'o  <uo 


(10) 
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where  the  subscript  o  denotes  conditions  at  the  free  surface.  Using  equations 
(9)  and  (10)  in  equation  (8),  we  obtain  the  expression  for  the  vorticity  at 
the  free  surface, 


Thus,  if  the  flow  in  the  interior  of  the  fluid  is  assumed  to  be  i rrotational , 
there  exists  a  boundary  layer  across  which  the  vorticity  reduces  to  zero. 
Note  that  such  a  boundary  layer  is  absent  if  the  free  surface  is  flat. 

If  surface  tension  is  neglected,  the  normal  stresses  at  the  free  surface 
must  be  constant  and  equal  to  the  ambient  pressure,  which  may  be  taken  to  be 
zero.  Equation  (6)  then  yields 


(IV) 

'9n'o 


=  0 


If  this  and  equation  (9)  are  used  in  the  continuity  equation  (1),  we  obtain 
the  rather  surprising  result 


i ,e.  u  =  constant  =  U 
o  00 

Since  the  velocity  at  the  free  surface  in  inviscid  flow  decreases  and  vanishes 

at  the  intersection  with  the  obstacle,  or  equivalently,  since  the  velocity 

outside  the  boundary  layer,  u.  say,  must  decrease  as  the  obstacle  is 

6 

approached,  the  above  result  is  physically  unrealistic  and  we  conclude  that 
the  condition  of  zero  or  constant  normal  stress  at  the  free  surface  cannot  be 
satisfied.  In  other  words,  the  influence  of  surface  tension  cannot  be 
ignored. 

Denoting  the  surface  tension  by  o,  the  balance  of  normal  stress  across 
the  free  surface  requires  that 

-  pa  +  ok  =  -  pQ  +  2  v  (  f£)0 


(12) 
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where  pa  is  the  (zero)  ambient  pressure  above  the  free  surface  and  p0  is  the 
pressure  in  the  liquid  just  at  the  free  surface.  Thus, 


( p  —  p  ) 

/9_V\  _  <  a  ra  Ko 

9n  'o  ”  2  p  ”  2  m 


(13) 


Substitution  of  equations  (9)  and  (13)  into  the  continuity  equation  (1)  now 
yields 


(—)  = 
^9s'o 


<  a  (Pa-Po* 

77  +  2-u 


(Pa’po) 


■jfs  (arctan  x)  +  - 2 


(14) 


where  is  the  slope  of  the  free  surface.  If,  as  a  first  approximation, 
we  assume  that  p0  =  pa,  i.e.  the  surface  tension  is  balanced  by  the  viscous 
stresses  rather  than  a  jump  in  pressure,  equation  (14)  can  be  integrated  to 
obtai n 


u„  =  U  -  tt"  arctan  r  (15) 

o  ®  2m  x 

The  distribution  of  vorticity  at  the  free  surface  is  then  given  by  equations 
(11)  and  (15)  as 


%=  "  2  K  V  '  2  K  (U»-  77  arctan 


(16) 


Equation  (15)  shows  that  the  velocity  along  the  free  surface  decreases  as 
increases.  The  previous  result  of  constant  velocity  in  the  absence  of 

surface  tension  is  recovered  when  a  =  0.  It  is  quite  surprising  that,  with 

the  above  approximation,  the  velocity  and  vorticity  at  the  free  surface  can  be 
predicted  solely  from  kinematics  (equations  1  and  9)  and  the  stress  conditions 
(equations  10  and  12)  at  the  surface  without  recourse  to  the  dynamical 

equations  (2)  and  (3).  Even  more  surprising  is  the  fact  that  equation  (15) 

embodies  a  SEPARATION  CRITERION  for  the  free-surface  boundary  layer,  since  uQ 
=  0  indicates  a  stagnation  point  on  the  free  surface.  Thus,  for  separation 


2uU 

=  tan  (-— ) 


sep 


*  ,  2W> 

=  tan  (-jjg) 


(17) 
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where  W  = 


and  Re  = 


are  the  Weber  and  Reynolds  numbers, 


a 


respectively,  based  on  some  characteristic  length  L.  Since  t,  and  depend 
upon  the  Froude  number,  F  =  U J  /gU  say,  equation  (17)  represents  a  separation 
criterion  in  terms  of  the  three  basic  nondimensional  parameters:  Vf,  Re  and  F. 


III.  SEPARATION  AT  THE  FREE  SURFACE  AND  BOW  VORTICES 


In  order  to  determine  if  equation  (17)  indeed  represents  a  plausible 
result,  we  seek  experimental  confirmation.  Consider  the  two-dimensional  flow 
ahead  of  a  semi -submerged  circular  cylinder  of  radius  a  with  its  axis 
horizontal  and  perpendi cul  ar  to  a  stream  of  velocity  U  in  the  positive  x- 


oo 


direction,  as  shown  in  figure  2.  Alternatively,  the  cylinder  moves  with 
velocity  U  in  the  negative  x-direction  in  a  liquid  at  rest.  The  latter 


oo 


corresponds  to  the  arrangement  in  the  experiments  of  Suzuki  (1975),  Honji 
(  1976)  and  Kayo  et  al .  (1982).  As  mentioned  in  the  Introduction,  a  vortex 
system  was  observed  ahead  of  the  cylinder  as  depicted  in  figure  2.  Tests  with 
different  velocities  indicated  a  wide  variation  in  the  length  6a  of  the 
vortex.  However,  the  three  sets  of  data  appear  to  be  in  some  conflict  with 
regard  to  the  influence  of  the  cylinder  speed  on  the  length  of  the  vortex. 
Although  Suzuki  appears  to  be  first  to  observe  the  vortices,  his  measurements 
are  related  to  the  breaking  of  bow  waves  ahead  of  the  cylinder  rather  than  the 
size  of  the  bow  vortex.  They  are  therefore  not  suitable  for  comparison  with 
the  present  theory.  The  measurements  of  Honji,  which  are  reproduced  in  figure 
3,  show  an  increase  in  6  with  increasing  Reynolds  and  Froude  numbers.  On  the 
other  hand,  the  observations  of  Kayo  et  al.,  shown  in  figure  4,  indicate  very 
large  and  scattered  values  of  6  at  the  lowest  velocities  in  the  tests.  With 
increasing  velocity,  6  appears  to  reach  a  minimum  and  show  a  moderate  increase 
thereafter.  Honji  also  showed  that  the  size  of  the  vortex  decreased  when  the 
surface  tension  was  reduced  by  adding  detergents  to  the  water.  Similar 
observations  were  also  made  by  Kayo  et  al.  but  no  quantitative  information  was 
obtained  for  the  case  of  the  horizontal  cylinder  considered  here. 


If  we  identify  the  most  upstream  point  S  of  the  vortex  (figure  2)  with 
the  free-surface  stagnation  point  predicted  by  the  present  theory,  then  it  is 
possible  to  calculate  the  distance  6a  from  equation  (17),  provided  the  slope 
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of  the  free  surface  is  known,  or  assumed,  ahead  of  the  cylinder.  The 
determination  of  the  free-surface  elevation  is  of  course  the  classical  problem 
of  nonlinear  ship-wave  theory. 

I II. 1  Evaluation  of  cx  from  inviscid,  irrotational-flow  theory 

Consider  a  circular  cylinder  of  radius  a,  half  immersed  in  a  uniform 
stream  of  velocity  U  in  the  positive  direction  of  the  x  axis.  Take  the 

oo 

origin  at  the  undisturbed  water  surface  and  the  y-axis  positive  upwards.  The 
irrotational-flow  velocity  components  in  the  x-  and  y-directions  are  U  and  V, 
respectively,  and  the  free-surface  elevation  due  to  the  presence  of  the 
cylinder  is  y  =  c.  (x) . 

The  condition  that  the  free  surface  is  a  streamline  gives 


(18) 


e  =  v/u 

x 


while  the  Bernoulli  equation  for  steady,  irrotational  flow  and  the  zero- 
pressure  condition  at  the  free  surface  yields  the  implicit  equation 


f(x,y)  =  1)2  +  V?  +  2gy  -  1) 2  =  0 


The  kinematic  boundary  condition,  Df/Dt  =  0,  then  yields  the  exact,  nonlinear 
boundary  condition 


U(UUX  +  VVX)  +  V(UUy  +  VVy  +  g)  =  0 


(19) 


where  the  subscripts  x,  y  denote  derivatives.  From  the  Cauchy-Riemann 
equations 


(20) 


equation  (19)  becomes 


V2UX  _  V(2UUy  +  g)  -  U2UX  =  0 


which  can  be  rewritten  in  the  nondimensional  form 
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V2UX  -  V( 2UUy  +  1/ ( 2F2) )  -  U2UX  =  0 


(21) 


where  (U,V)  and  (x,y)  are  scaled  with  U  and  a,  respectively,  and  the  Froude 

oo 

number  F  is  defined  by 


U 

oo 

/2ga 


(22) 


Combining  equations  (18)  and  (21),  we  obtain 


-  4F2UUx 

Cx  = - 5 - 2 - 2 - 2  1/2 - 2 —  (23) 

[(1+4F  UU)  +  (4F  UUj  ]  '  +  1  +  4F  UU 

y  x  y 

which  represents  the  free-surface  slope  upstream  of  the  cylinder.  Since, 
however,  the  velocity  field  and  its  derivatives  are  difficult  to  determine, 
especially  close  to  the  body,  their  values  from  the  double-body  approximation, 
i  .e., 


U 


1 


2  2 
x  +  y 


+ 


il 


/  2  2.2 

(x  +  y  ) 


V  = 


_ 2  x  y 

.  2  2,2 

(x  +  y  ) 


(24) 


evaluated  at  y  8  0,  are  used.  The  reason  for  using  the  double-body  solution 
rather  than  that  from  linear  wave  theory  is  that  the  former,  which  considers 
the  effect  of  nonlinearity,  gives  a  better  approximation  to  the  velocity  field 
near  the  body.  A  graph  of  as  a  function  of  x/a  is  given  in  figure  5  for 
several  values  of  the  Froude  number. 


Substitution  of  £x  from  equation  (23)  with  U,  Ux  and  Uy  at  y  =  0  from 
equation  (24)  into  equation  (17)  gives 

8F2[(1  +  e)2-  1]  _  t  2W  ,oq, 

{(1  +  3)  +  64F  [(1  +  3)  -  1]}  +  (1  +  3) 

which  can  be  solved  numerically  (or  graphically  using  figure  5)  for  the 
location  of  the  separation  point  3a  (defined  in  figure  2)  for  given  values  of 
F,  Re  and  W.  When  the  equation  has  several  positive  roots,  the  largest  one  is 
taken  since  separation,  if  present,  would  occur  at  the  most  upstream  point. 
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111.2  Comparison  with  experimental  results 

Equation  (25)  has  been  solved  for  the  conditions  in  the  experiments  of 
Honji  (a  =  0.05  m)  and  Kayo  et  al.  (a  =  0.10m).  For  Honji's  data,  the  value 
of  y  (=  1.22  x  10"3  Ns/m?)  was  inferred  from  the  quoted  Reynolds  numbers  and 
velocities.  For  the  data  of  Kayo  et  al.,  their  quoted  value  of  0.913  x  10" 3 
Ns/m?  was  used.  In  both  cases,  a  is  assumed  to  take  the  standard  value  of 
0.07  3  N/m.  The  results  are  shown  by  the  dotted  lines  in  figures  3  and  4. 
Note  that,  in  both  figures,  the  Reynolds  and  Froude  numbers  have  been 
redefined  consistently  on  the  basis  of  the  cylinder  diameter,  and  the  scales 
adjusted  accordingly. 

The  results  for  Honji's  experiments,  which  were  conducted  at  very  low 
Reynolds  and  Froude  numbers,  indicate  qualitative  agreement  between  theory  and 
observations,  but  the  predicted  delayed  separation  indicates  smaller 
vortices.  In  the  case  of  Kayo  et  al.,  the  theory  predicts  the  trends  at  the 
higher  Reynolds  and  Froude  numbers  but  the  calculations  now  indicate  earlier 
separation  and  larger  vortices.  The  experimental  results  at  very  low  Froude 
numbers  show  considerable  scatter.  At  the  lowest  Froude  number  in  these 
tests,  the  experiments  indicated  very  large  separation  distances  which  varied 
from  test  to  test  over  a  wide  range.  The  authors  noted  that  the  largest 
distances  were  observed  during  the  first  experiment  in  the  morning!  The 
theory  obviously  does  not  predict  this  phenomenon. 

In  summary,  we  note  that  the  present  simple  theory  appears  to  provide  an 

explanation  for  the  origin  of  the  bow  vortices.  The  contradictions  in  the 

available  experimental  data  indicate  a  need  for  more  refined  and  controlled 

experiments.  On  the  other  hand,  the  disagreement  between  the  theory  and 

existing  data  could  be  due  to  a  number  of  factors  which  need  to  explored 

further.  Among  these  are  (a)  the  influence  of  the  assumption  of  no  pressure 

jump  across  the  free  surface,  (b)  the  uncertainty  in  the  determination  of  the 

free-surface  slope  t,  ,  (c)  the  interaction  between  the  bow  vortex  and  the 

x 

shape  of  the  free  surface,  (d)  the  occurrence  of  turbulent  flow  ahead  of 
separation,  and  (e)  the  influence  of  surface  contaminants  in  the  experiments. 
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IV.  FREE- SURFACE  BOUNDARY-LAYER  EQUATIONS 

The  results  presented  above  have  been  obtained  solely  from  an  examination 
of  the  boundary  conditions  at  the  free  surface  and  the  assumption  that  there 
is  no  pressure  jump  across  the  free  surface.  In  order  to  remove  the  latter 
restriction  it  is  necessary  to  seek  a  solution  of  the  equations  of  motion. 

The  Navier- Stokes  equations  may  be  simplified  to  obtain  the  equations  of 
the  free-surface  boundary  layer.  For  this  purpose,  it  is  necessary  to  make  an 
order-of-magnitude  analysis  utilizing  the  known  boundary  conditions.  Let  L  be 
a  characteristic  length  and  6  be  the  boundary-1 ayer  thickness.  Then  K*  =  <L, 
the  nondimensional  free-surface  curvature  which  gives  rise  to  the  boundary 
layer,  is  controlled  by  the  body  geometry  and  the  Froude  number. 

The  equation  of  continuity  (1)  and  its  derivative  are 


3u  av 

3s  +  hl  an  +  <v  =  0 

2  2 

a_ju_  ,  a__v  0  a_v  _  n 

asan  +  hi  2  +  2k  an  ~  0 
3  n 


Since  hj  =  0(1),  from  the  former,  we  have 


3  v 


y-n  -  0  O  -  0  0)  IT 


(2b) 

(27) 


(28) 


and  then 


v  -  0  (ft  U„ 


Also,  equation  (10)  shows  that 
U 


yw=  o  (K  ) 


"T 


which  gives 


2  *  II 

3  U  _  n  rK  L1  u- 

~2  ~  0  (  T-'  r? 

an  L 


(29) 


(30) 


(31) 


From  equations  (27)  and  (30),  we  obtain 
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9-4  =  0(K*) 
3  n 


U 

00 


(32) 


Now,  since  the  first  convective  term,  together  with  the  pressure  and 
gravity  terms  in  the  momentum  equation  (2),  yields  the  Bernoulli  equation,  we 
conclude  that  the  leading  viscous  term  must  be  of  order  (<6),  i.e. 


v 


32u 

77 


=  0 


(U) 


0(<6 ) 


oo 


r~ 


Hence 


f  * 0  (-b 

L  /Re 

and  equations  (29)  and  (31)  become 


(33) 


v  *  0  u. 

/Re 

^2  _  U 

— *  =  0  (KVRe)  ~ 

3n  L 

Thus,  we  have  the  following  results: 


(34) 

(35) 


u  -  0  (1) 

9u  _  n  ,  *»  U“ 

an  -  0  (k  )  r~ 

2  U 

i-lj  =  0  (KVRe) 
3n^ 


7 


*  =  0  (V  u 

/-K -  00 

/l?e 

*  U 

I"  =  0  (K  )  ” 

3n  L 


(36) 


When  the  above  orders  of  magnitude  are  introduced  in  the  Navier-Stokes 
equations  (1)  -  (3)  and  terms  up  to  order  (k<5)  are  retained  in  each  equation, 
we  obtain  the  boundary-layer  equations 


iH  +  i*  .  o 

3S  8tl  u 


(37) 


32u 

”  T? 

3n 


(38) 
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-  ku2  +  fn  (p  +  9y)  =  0  (  39) 

IV. 1  Velocity  Profile 

The  known  boundary  conditions  can  also  be  used  to  determine  some 
characteri sties  of  the  velocity  distribution  across  the  boundary  layer.  To 
the  first  order,  the  vorticity,  from  equation  (8),  is 

“  f  -  TF  "  KU  (40) 

The  condition  of  zero  vorticity  outside  the  boundary  layer  yields 


(— )  = 
v3n;S  • 


K  Ut 


(41) 


For  positive  < ,  equations  (10)  and  (41)  show  that  the  velocity  reaches  a 
maximum  within  the  boundary  layer. 

A  velocity  profile  which  explicitly  satisfies  the  boundary  conditions 
(10)  and  (41)  is  the  cubic 

rii  .  teex  2  .  ..  3~i  f  a  o\ 

u  =  u0[l  +  <n  +  n  +  — ^  n  ]  (42) 

6 


where 


o  =  -(l+a)(2+ic5) 


_  ^  ( 1  +  a ) 


(43) 


and  a  is  a  parameter  which  ensures  compatibility  with  equation  (2)  at  the  free 
surface,  i.e. 

r\  I 

6  ,3  ux  6  ,  ,  i  .  Ox 

a  =  W  ^  2^0  =  2^iT  Kuo  +  9C  + 

o  3  n  o 

where  the  primes  denote  derivatives  with  respect  to  s.  The  pressure  p0  can  be 
eliminated  from  the  above  using  the  normal -stress  condition  (12)  and  the 
equation  of  continuity  (37)  since 


p1  =  -  a< 1 

Ko 


(44) 
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Thus, 


6 

2vku 


KU0 


gc  - 


OK 


~  2vu") 


(45) 


It  should  be  noted  that  the  combination  (uQu0'  +  gc ' )  is  small  and  therefore 
a  is  of  order  unity. 

From  equation  (42),  the  ratio  of  velocity  at  the  edge  of  boundary 
layer,  ufi ,  to  that  at  the  free  surface,  u  ,  can  be  expressed  in  terms 
of  a  as 


\  =  3  +  <5 ( 2+ct )  f 

U  3  +  K  6  • 

o 

This  indicates  that  a  > 


1  +  (l+a) 

-  1  for  u„  >  u  . 

5  o 


(46) 


IV. 2  Momentum  Integral  Equations 

As  a  first  step  in  the  solution  of  the  first-order  boundary-layer 
equations  ( 37 ) - ( 39 ) ,  it  is  convenient  to  use  the  well-known  integral  method. 
Integration  of  the  normal  momentum  equation  (39)  yields 

(p  +  gy)  -  (jp  +  gc)  +  /  <u2dn 

0 

and  therefore 

(f  +  gy)6  =  (jp  +  gc)  +  /  <u2dn 

and 

(^  +  gy)  =  -gi-  (^  +  gc)  +  Is  I  Ku2dn  (49) 

0 

Now  applying  the  Bernoulli  equation  at  the  edge  of  the  boundary  layer,  we  have 


(47) 


(48) 


1 

2 


+  (^  + 


1  2 
9 y)6=  2  u» 


(50) 


Substitution  of  this  in  equation  (48)  gives 
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(~  +  94)  =  7  (U~  '  uf)  -  J  <u2dn 
and  from  equation  (49) 

Is  (f  +  gy)  =  -  u6uj  -  {/  <u2dn} 

n 

Also,  integration  of  the  continuity  equation  (37)  gives 


(51) 


(52) 


v 


n 


/ 

0 


au 

Ts 


dn 


(53) 


Using  equations  (52)  and  (53),  together  with  the  boundary  conditions  of 
equations  (10)  and  (41)  in  the  viscous  term,  equation  (38)  can  be  integrated 
across  the  boundary  layer  from  n  =  0  to  n  =  5.  After  some  rearrangement  this 
yields  the  momentum-integral  equation  in  the  form 


ds  u(u  -  u«)d"  *  u«  SQ  (u  -  V  dn 

6  n  d  6  6  2 

-  /  <u  vr  /  udndn  -  -rr  /  /  <u  dndn  +  v<(urt  +  uj 

0  35  0  ds  0  n  06 


=  0 


(54) 


If  the  integrals  in  the  above  are  evaluated  using  the  velocity  profile  of 
equation  (42),  and  terms  of  0(c6 )  are  neglected,  we  obtain  the  rather  simple 
momentum  integral  equation 


(55) 


so  that 

62  = - “—7-  /  <uods  (56) 

(  3  +  a  )  UQ<  -<*> 

Since  a  is  given  by  equation  (45),  the  above  equation  relates  the  boundary- 
layer  thickness  to  the  free-surface  velocity,  uQ,  and  shape,  r,  and  <. 

Another  equation  containing  the  above  variables  can  be  obtained  from  the 
normal -momentum  equation  (51)  by  substituting  into  it  the  velocity-profile 
relations  (42)  and  (46),  and  using  equations  (44)  and  (45).  This  leads  to 
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6  lls  {3  (4  +  a)  uo  <6}  =  '  (57) 

Thus,  equations  (56)  and  (57)  can  be  solved  for  g  and  uQ  if  the  shape  of 
the  free  surface  (i.e.,  c  and  <)  is  prescribed.  Note  that  surface  tension  is 
involved  through  the  parameter  a,  and  separation  would  be  indicated  by  uQ  =  o. 

It  is  also  of  interest  to  evaluate  the  displacement  and  momentum 
thicknesses  of  the  boundary  layer.  If  the  flow  were  potential,  the  velocity 
would  vary  according  to  equation  (40)  with  w  =  0.  This  gives  the  potential- 
flow  velocity  (up)  distribution 

up  =  u6{ 1+k (6-n)}  (58) 

The  integral  thicknesses  defined  by 
*  1  6 

<5  =  —  /  (u  -  u)dn  (59) 

6  o  H 

and 

0  =  — ^  f  u(u  -  u)dn  (6°) 

u6  o 

can  then  be  evaluated  using  equations  (42),  (46)  and  (58)  to  obtain 

<6  =  -g-  (3  +  a)  (k6)^  +  0(k6)^  (61) 

and 

<0  =  ■g'  (3  +  a)  (<6)^  +  0(<6)^  (62) 

If  equation  (62)  is  combined  with  equation  (56),  we  obtain  the  rather  simple 
resul t 

e  =  -^  /  <u2ds  (63) 

u3  0 

0 

Since  equation  (63)  does  not  contain  the  parameter  a,  it  can  be  readily 
integrated  with  the  condition  0  =  0  at  s  =  -  <*>,  provided  <  and  u0  are  known. 


i 
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to  obtain  an  estimate  of  the  momentum  thickness.  As  a  first  approximation,  we 
use  the  inviscid,  i rrotational -fl  ow  analysis  of  Section  III,  i.e.,  equation 

(23)  to  determine  <  and  equation  (24)  to  obtain  u0  (=  y  =  0),  Figure  6 

shows  the  development  of  the  momentum  thickness  at  Re  =  2  x  10^  for  several 
Froude  numbers.  It  is  seen  that  9/a  is  of  the  order  of  10"^  for  distances  of 
the  order  of  one  diameter  upstream  of  the  cylinder.  The  rapid  growth 
of  9  ahead  of  the  obstacle  also  indicates  the  increased  proneness  to 
separation. 


V.  CONCLUDING  REMARKS 

The  flow  ahead  of  a  semi -submerged  two-dimensional  body  moving  through  a 
viscous  liquid  is  considered.  Examination  of  the  exact  boundary  conditions  at 
the  free  surface  has  shown  that  surface  tension  plays  a  critical  role  in  the 
determination  of  the  free-surface  velocity.  The  simplifying  assumption  of  no 
pressure  jump  across  the  free  surface  then  leads  to  the  prediction  of  a 
stagnation  point  ahead  of  the  body.  This,  in  turn,  explains  the  origin  of 
vortices  ahead  of  the  body  which  have  been  observed  in  experiments.  the 
predictions  of  the  simple  theory  are  in  qualitative  agreement  with 
experimental  results.  Possible  reasons  for  some  discrepancies  in  the 
available  experimental  data,  and  the  lack  of  more  precise  agreement  with  the 
theory,  have  been  suggested. 

A  detailed  analysis  of  the  equations  of  motion  has  been  carried  out  to 
derive  the  equations  of  the  free-surface  boundary  layer  which  contains  the 
vorticity  generated  by  the  surface  curvature.  An  integral  method  has  been 
utilized  to  obtain  two  ordinary  differential  equations  which  relate  the 
boundary-layer  thickness  and  free-surface  velocity  to  the  curvature  and  slope 
of  the  free  surface.  These  remove  the  restrictive  assumption  on  the  pressure 
jump  and  should  lead  to  a  more  accurate  prediction  of  the  separation  and 
vortex  location.  Solutions  of  these  equations  are  in  progress. 

Finally,  we  note  that  the  analysis  presented  here  is  restricted  to  a  two- 
dimensional  flow.  Nevertheless,  it  provides  an  explanation  for  the  existence 
of  vortices  observed  ahead  of  ship  models.  Also,  it  can  be  generalized  for 
application  to  three-dimensional  flows  and,  in  combination  with  the  kinematic 
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theories  of  vorticity  amplification  (Hawthorne,  1954;  Lighthill,  1956),  as 
suggested  by  Mori  (1984)  and  Takekuma  and  Eggers  (1984),  may  lead  to  a 
rational  theory  for  the  prediction  of  the  necklace  vortices  ahead  of  ships. 
The  complex  interaction  between  the  free-surface  boundary  layer,  the  bow 
vortices  and  the  overall  flow  pattern  around  the  bow  is  a  subject  for  further 
research,  as  is  the  connection  between  the  bow  vortices  and  the  breaking  of 
bow  waves. 
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FIGURE  1.  NOTATION 


FIGURE  2.  VORTICES  AHEAD  OF  A  TWO-DIMENSIONAL  CYLINDER 
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FIGURE  3.  EXPERIMENTS  OF  HONJI  (1976) 
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FIGURE  4. 


EXPERIMENTS  OF  KAYO  ET  AL.  (1982) 
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FIGURE  5.  FREE-SURFACE  SLOPE  AHEAD  OF  A  CYLINDER 
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FIGURE  6.  MOMENTUM  THICKNESS  OF  THE  FREE-SURFACE  BOUNDARY  LAYER  AHEAD  OF 
THE  CYLINDER 
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